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ABSTRACT 

We report a significant detection of weak, tangential distortion of the images of cosmo- 
logically distant, faint galaxies due to gravitational lensing by foreground galaxies. A mean 
image polarisation of (p) = 0.011 ± 0.006 (95% confidence bounds) is measured for 3202 
pairs of source galaxies with magnitudes in the range 23 < r < 24 and lens galaxies with 
magnitudes 20 < r < 23. The signal remains strong for lens-source separations <J 90 arc- 
sec, consistent with quasi-isothermal galaxy halos extending to large radii (<; 100/i -1 kpc). 
Our observations thus provide the first evidence from weak gravitational lensing of large 
scale dark halos associated with individual galaxies. The observed polarisation is also 
consistent with the signal expected on the basis of simulations incorporating measured 
properties of local galaxies and modest extrapolations of the observed redshift distribu- 
tion of faint galaxies. From the simulations we derive a best-fit halo circular velocity of 
V* ~ 220 km s _1 and characteristic radial extent of s* <; 100 h~ x kpc. Our best-fit halo 
parameters imply typical masses for the lens galaxies within a radius of 100 h~ x kpc on the 
order of l.OiJ'y X 10 12 h~ x Mq (90% confidence bounds), in good agreement with recent 
dynamical estimates of the masses of local spiral galaxies. This is particularly encouraging 
as the lensing and dynamical mass estimators rely on different sets of assumptions. Con- 
tamination of the gravitational lensing signal by a population of tidally distorted satellite 
galaxies can be ruled out with reasonable confidence. The prospects for corroborating and 
improving this measurement seem good, especially using deep HST archival data. 
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1. INTRODUCTION 



The notion that cosmologically-distributed masses might cause weak but measurable 
changes to the shapes of distant galaxies has a long history. The most striking example 
of this "cosmological distortion effect" (Kristian & Sachs 1966) is the distortion of distant 
galaxies into giant arcs caused by rich clusters of galaxies (eg. Fort & Mellier 1994). At- 
tempts have also been made to measure the weak distortion of distant galaxy images by 
mass fluctuations associated with large-scale structure (eg. Kristian 1967; Valdes, Tyson 
& Jarvis 1983; Mould et al. 1994, Paper I). These studies have progressed to the point 
where they can produce interesting limits on the large scale distribution of mass in the 
universe. The purpose of this paper, however, is to show how weak gravitational lensing 
can be used to study the mass distributions on much smaller scales, those associated with 
individual galaxies. 

Galaxies constitute the fundamental building blocks of luminous structure in the uni- 
verse, yet we are largely ignorant of such basic physical parameters as their typical masses 
and radial extents. Popular theories of galaxy formation predict that most bright galax- 
ies should reside in massive (<; 1O 12 M0) dark halos that extend far beyond their optical 
radii (<; 100/i -1 kpc). Direct observational evidence to test such theories is, however, 
scarce. While observations of the central parts of galaxies provide good mass estimates 
for these regions, the lack of information on the form and extent of the dark halos of 
individual galaxies limits our determination their masses. Observations of local galax- 
ies, most importantly our own (eg. Fich & Tremaine 1991), indicate that the majority of 
bright spiral galaxies have dark halos which extend isothermally out to at least ~ 30 kpc. 
Studies using the dynamical properties of ensembles of faint companions to samples of 
bright spirals favour a continuation of the halo out to radii ~ 100 kpc (Zaritsky & White 
1994). On purely theoretical grounds, if we suppose the mass of the halos increases linearly 
with radius out until their density reaches the critical density, they will have outer radii 
~ (2 /3) 1 1 2 (V c I Ho) ~ 1 — 2/i -1 Mpc, where V c is the circular velocity. Alternatively, the 
halos of most giant galaxies may be truncated at radii smaller than a tenth of this value. 

In this paper we will report on a significant detection of distortion of the shapes of 
distant galaxy images due to weak gravitational lensing by individual foreground galaxies. 
Using the observed gravitational lensing signal we will then attempt to measure masses 
for typical field galaxies. The advantages of the lensing approach to determining galaxy 
masses are two-fold. Firstly, it is capable of probing the halos of galaxies out to very large 
radii, r <; 100/i -1 kpc, where few classical techniques are viable. Secondly, the lensing 
analysis is relatively unaffected by the dynamical properties of the possibly unvirialised 
outer regions of the halos. Lack of detailed information about the dynamics at large radii 
in galaxies is one of the central problems for the application of dynamical mass estimators. 
The dynamical and lensing mass estimates depend upon different model assumptions and 
a comparison of the results obtained from the two techniques is one of the few ways in 
which the validity of these assumptions can be tested. 

The signal for which we are searching is a distortion of the images of faint galaxies 
resulting in a weakly preferred alignment of faint galaxies around brighter galaxies. If the 
faint galaxies are gravitationally lensed by the brighter systems, the major axes of their 
images will tend to lie perpendicular to the radius vectors joining the faint galaxies to the 
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Fig. 1: Orientation of faint galaxies relative to bright galaxies 

lens centres (Fig. 1). The strength of this signal depends upon the distances of the lens 
and source galaxies, the mass of the lens, and the angular separation of the lens and source 
on the sky. With statistical information from direct spectroscopic studies (eg. Lilly 1993, 
Tresse et al. 1993) on the redshift distributions of both the source and lens galaxy samples 
we can therefore solve for the typical masses of the lensing galaxies. 

The detection of galaxy-galaxy lensing was first attempted by Tyson et al. (1984) us- 
ing scans of photographic plates (cf. also Webster 1983) from which they tried to measure 
an excess of background galaxy images tangentially elongated with respect to brighter, can- 
didate lens galaxies. The galaxy sample available to Tyson et al. was very large, ~ 47, 000 
background galaxies with 22.5 < J < 23.5 and ~ 12,000 lens galaxies with 19 < J < 21.5. 
Faint image position angles were measured for over 28, 000 galaxy pairs and no statisti- 
cally significant difference from an isotropic distribution was seen for galaxy separations 
greater than ~ 3 arcsec. On this basis, Tyson et al. concluded that the typical galaxy 
circular velocity had a surprisingly small upper bound, <J 170 km s _1 . Unfortunately, the 
relatively modest seeing of their plate material could have a very detrimental effect on 
the efficiency of the test. Moreover, Kovner & Milgrom (1987) performed a more careful 
calculation of the magnitude of the effect anticipated, taking into account integration over 
galaxy luminosity functions and distances as well as the correlation of internal, galaxian 
velocity dispersions with luminosity and concluded that the observation of Tyson et al. 
was consistent with conventional dynamical models of local galaxies. 

In §2 we describe our observations and present our results for the detection of weak 
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gravitational lensing of distant galaxy images by foreground galaxies, together with the 
analyses we have performed to attempt to understand the importance of systematic errors. 
In §3 we describe model calculations, both analytic and Monte Carlo, that we have used to 
translate our measured galaxy-galaxy lensing signal into a quantitative statement about 
galaxy masses and extents, results of which are presented in §4. Finally, in §5 we relate 
our measurement of galaxy masses and radial extents to more conventional determinations 
of these quantities. 

2. MEASUREMENT OF GALAXY-GALAXY LENSING 

To estimate roughly the expected strength of the galaxy-galaxy gravitational lensing 
signal, let us model a lens galaxy as a singular isothermal sphere with circular velocity V c . 
An ellipticity ~ 27rV c 2 /c 2 6 will then be induced in the image of a background faint galaxy 
located an angular distance ~ 6 from the lens. This is of the order of a few percent effect 
for galaxy pairs with separations 6 ~ 30 arcsec where the lens galaxy is a typical bright 
spiral galaxy. In order to detect this signal in the presence of the noise associated with 
the intrinsic galaxy shapes, over a thousand foreground-background galaxy pairs must be 
measured. If a sufficiently large number of pairs are available, it may also be possible to use 
the dependency of the lensing signal on the distance from the lens centre, to study the 
angular extent of galaxy halos. In the following section we discuss an observational dataset 
which should be of sufficiently quality, depth and size to allow us to detect galaxy-galaxy 
lensing. 

2.1 Observational Data 

The imaging data used in our analysis is of a single 9.6 X 9.6 arcminute blank field 
centred on a(1950) = 17 ft 21 m 07 s (5(1950) = +49°52'21", taken in Gunn r. The data 
were acquired in periods of good seeing, 0.7-0.9 arcsec, using the direct imaging mode 
of the COSMIC imaging spectrograph (Dressier et al. 1995) on the 5-m Hale telescope, 
Palomar. These data have been used previously to study the coherent distortion of faint 
galaxy images due to weak gravitational lensing by large-scale structure (Paper I) and 
the angular clustering statistics of faint galaxies (Brainerd, Smail, & Mould 1995). The 
reduction of the data to a catalogue of detected objects is detailed in Paper I. 

The final stacked r image used for our principal data analysis consists of a total of 
24.0 ksec integration, has a la surface brightness limit of /i r = 28.8 mag arcsec -2 , seeing 
of 0.87 arcsec FWHM, and a total area of 90.1 arcmin -2 . The object catalogue created 
from this frame using the FOCAS image analysis package (Valdes 1982) contains ~ 6600 
objects brighter than the 80% completeness limit of r = 26.2. Adopting a conservative 
magnitude limit of r < 26.0, where the detections are approximately 97% complete, we 
obtain a cumulative surface density of 71.8 galaxies arcmin -2 or 2.6 X 10 5 deg -2 . Due 
to the presence of classical distortion in the corners of the frame, our analysis uses only 
those galaxies which lie within a circle of radius 4.8', centred on the chip. There are 4819 
galaxies within this area brighter than r = 26.0. 

An additional, shallower (/-band image of the same field is also available and is used 
solely for the purpose of providing colour information on the objects detected in the deep 
r image. The (/-band image has a total exposure time of 6.0 ksec and seeing of 1.2 arcsec. 
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The median colour errors for galaxies detected on the r image with (g — r) < 1 at r ~ 26 
(approximately 80% of the total population) are A(g — r) ~ 0.2. 

In order to calculate the gravitational lensing signal yielded by our models and to 
estimate the lensing-induced mean image polarisation of our faint objects, we shall require 
the distribution function of the intrinsic source galaxy ellipticities. To linear order the 
source ellipticity distribution can be estimated by the ellipticity distribution of the images. 
We have, therefore, measured the distribution function of image ellipticities in our sample 
and find that it is adequately fit by the normalized distribution function 

P e (e) = 64eexp[-8e] (2.1) 

with mean ellipticity (e) = 0.12 and we estimate the error to be ±0.02. 
2.2 Position angle probability distribution 

In this section we investigate the orientation of faint galaxies relative to the directions 
to nearby bright galaxies (see Fig. 1). The unweighted second moments of the intensity 
were measured and a complex orientation, x, formed for each faint image. The modulus 
of x is (° 2 ~~ b 2 )/(a 2 + b 2 ) where a 2 and b 2 are the principal second moments of the 
intensity provided by FOCAS. (See Paper 1 for a discussion of alternative prescriptions 
for describing the image shapes.) For a ~ 6, \x\ is approximately equal to the measured 
ellipticity, e = 1 — 6/a, and this identification is adequate for our purposes. The phase of 
X is twice the position angle, </>, and we shall use the convention of measuring <f> to be the 
angle between the major axis of the equivalent ellipse and the radius vector, measured in 
a counter-clockwise sense. We also combine positive and negative position angles so that 
<f> is restricted to [0,7r/2]. 

In the absence of distortion of the faint galaxy images, we expect a uniform distribution 
of their position angles for all projected separations, 6, between the candidate lenses and 
sources. When the background galaxy images are gravitationally lensed by foreground 
galaxies, we expect a distribution of position angles that is non-uniform with a deficit of 
faint images oriented radially ((f) = 0) and an excess of faint images oriented tangentially 
((f) = vr/2). As we show in §3, the deviation of the distribution from uniform should exhibit 
a cos 2(f> variation. At large projected separations, of course, the distribution of (f> should 
become uniform. 

In Fig. 2 we show the observed faint galaxy orientation distribution, P^(^), for lens- 
source projected separations of 5 < 6 < 34 arcsec. To determine P<f>((f>) we use the 439 
bright galaxies with magnitudes in the range 20.0 < r& < 23.0 (the candidate lenses) and 
faint galaxies with magnitudes in the ranges: (a) 23.0 < r s < 24.0, (b) 23.0 < r s < 25.0 
and (c) 23.0 < r s < 26.0. The magnitude limits used to select the candidate lens galaxies 
are roughly equivalent to the depth of current faint spectroscopic surveys (eg. Lilly 1993, 
Tresse et al. 1993), which provide statistical distance information on this population. The 
number of faint images and faint -bright pairs as a function of limiting source magnitude 
are summarised in Table 1. The lower limit on 6 was chosen so as to avoid overlapping 
faint and bright image isophotes. We discuss the choice of upper limit on 6 below; the 
value here corresponds to an average impact parameter of ~ 120/i -1 kpc at our median 
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Fig. 2: Probability distribution P(f>(<f>) of orientation of faint galaxies relative to the direc- 
tions of bright galaxies with projected separations 5 < 6 < 34 arcsec. The bright galaxies 
have 20 < r < 23 and the faint galaxies have magnitudes in the ranges indicated in the 
text and figure panels. For the best case, (a), of a non-uniform P^,(</>), the best fitting 
theoretical cos 2<f> variation is also shown. 



lens redshift. The error bars on P(f>(<f>) are estimated from bootstrap resampling of the 
data. 

From Fig. 2, the (binned) distribution of faint image position angles appears non- 
uniform for the two brighter galaxy samples. In addition, P(f>(<f>) for the two brighter 
samples is qualitatively in agreement with what one would expect in the case of lensing by 
the bright foreground galaxies. We investigate the significance of the deviation of P(f>(<f>) 
from a uniform distribution using standard statistical techniques. 

A x 2 test performed on the binned P(f>(4>) in Fig. 2 rejects a uniform distribution 



6 



(P^(^) = 2/tt) at a modest confidence level (> 97%) for both of the brighter samples 
and cannot reject a uniform distribution for the faintest sample. The sensitivity of the 
X 2 test to the arbitrary choice of bin size when using a continuous variable is well known. 
With this in mind, we consider the Kolmogorov-Smirnov (KS) test, which compares the 
observed continuous, cumulative distribution with a similar theoretical distribution, to be a 
better statistic for our purposes. KS tests performed on the cumulative P(f>(<f>) distributions 
corresponding to the binned distributions shown in Fig. 2 reject uniform distributions at 
better than the 99% confidence level for both brighter source samples, however we still can 
not rule out the hypothesis that our faintest sample is uniformly distributed. A summary 
of rejection confidence levels as a function of the source limiting magnitude is given in 
Table 1. 

We do not view the absence of an obvious lensing signal in our faintest source sample 
as a strong concern. In a universe where galaxies at fainter magnitudes are observed 
at progressively more distant epochs we might expect the lensing signal to continue to 
strengthen in fainter samples. Unfortunately, as we reach fainter in our Universe the 
typical angular size of galaxies decreases, and by r ~ 26 they have half-light radii of <> 0.4 
arcsec. This is small compared to our seeing and means that we only partially resolve our 
faintest galaxies, leading to increasing errors in the determination of their ellipticities and 
position angles, degrading the gravitational lensing signal that we measure. Such a trend 
is clear from the data exhibited in Fig. 2. We must, therefore, settle for an optimal choice 
of limiting magnitude for the faint images and, from Fig. 2, this appears to be r ~ 24 for 
our data. We therefore adopt a source sample defined by 23.0 < r s < 24.0 (Fig. 2a) for 
our principal analysis. 

As we anticipate P(f>(<f>) will exhibit a cos 2<f> variation in the case that the faint images 
have indeed been gravitationally lensed by the brighter galaxies, we show in Fig. 2a the 
best fitting cos 2<f> variation for P(f>(<f>) for the subsample that will be used for our principal 
analysis. 

2.3 Possible systematic effects 

We have performed a number of tests to investigate possible systematic effects in 
our data which would give rise to the observed non-uniform P^cf)). To begin with we 
considered lens-source pairs with projected separations of 5 < 6 < 34 arcsec, lens galaxies 
with 20 < rd < 23, and source galaxies with 23 < r s < 24, and computed P(f>(<f>) for the 
following: 

(1) <f> taken to be the position angle of the lens galaxies relative to the lines connecting 
their centroids with those of the sources 

(2) the position angle of the source images relative to lines connecting their centroids with 
random points (i.e. not corresponding to the centroids of the lenses) 

(3) the position angle of the source images relative to the lines connecting their centroids 
to bright stars on the frame 

(4) random position angles were substituted for the observed position angles of the faint 
images 

(5) using the observed centroids and source position angles, the distributions of positive 
and negative values of <f> were computed independently and compared 
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In cases (l)-(4), neither the \ 2 nor the KS test rejects the uniform distribution (i.e. no 
"signal" is observed in these cases). In case (5), the distributions obtained for the two 
subsamples are statistically indistinguishable and the KS test rejects uniform distributions 
at the ~ 99% confidence level (a somewhat weaker rejection than that obtained when 
positive and negative values of cf> were combined, cf. Table 1). 

Next we discuss the possible role of a point spread function (psf ) asymmetry in creating 
the observed signal. The presence of a non-circular psf or guiding errors introduces a 
preferred orientation in the object images and can give rise to a non-uniform P(f>(<f>) on 
very large angular scales. When the orientation of the faint galaxies relative to the bright 
galaxies is computed on small scales (<^ 1/4 the size of the frame), this effect is canceled 
out due to the fact that we can compute the orientations of the faint images in complete 
annuli around the majority of the bright centers. However, on scales ^ 1/2 the size of the 
frame, the effect becomes very significant as we can no longer compute the orientations of 
the faint images in complete annuli around most of the centers. In general, the edge effect 
at large scales that gives rise to the non-uniform P(f>(<f>) will not resemble that expected 
for a gravitational lensing signal, being of a larger amplitude and peaking at an angle cf> 
corresponding to a combination of the preferred orientation of the faint images and their 
direction vectors relative to the bright centers. 

In the case of our data, there is a slightly elliptical psf measured from the bright stars 
on the frame and a corresponding weakly preferred image orientation. Computing P(f>(<f>) 
for our fiducial subsamples of galaxies, we find it is consistent with a uniform distribution 
on scales of ~ 100-150 arcsec (as would be expected for a gravitational lensing signal), 
but on scales <; 250 arcsec, P(f>(<f>) is significantly non-uniform due to the edge effect. The 
form of P(f>(<f>) on these scales, unlike Fig. 2, does not coincide with the expectations of 
a gravitational lensing signal at these scales, being of a larger amplitude and peaking at 
cf> = 7r/4 with corresponding suppressions at cf> = and cf> = tt/2. Such a signal (both 
amplitude and shape) in our data is expected at this scale given the observed preferred 
orientations of the objects (mean position angle of ~ — 57r/18). 

Finally, as a check for systematics associated with the image detection routines and 
second moment determinations, P(f>(<f>) was computed from a simulated data frame in which 
a signature of weak gravitational lensing was not included. Images on the simulated frame 
were assigned random locations and orientations. The galaxy parameters (ellipticities, scale 
sizes, and magnitudes) were chosen such that after the images were convolved with the 
telescope psf, the distribution of these parameters matched the observed distributions. The 
surface density of objects above the completeness limit was matched to that observed and, 
in addition, fainter galaxies were added by extrapolating the number counts 2 magnitudes 
fainter in order to simulate the effects of crowding and merging of these undetected faint 
galaxies on the detected objects. The pixellated images were convolved with the psf defined 
by the telescope and atmosphere, and sky noise was added to obtain the same detection 
limits as in the observed images. The simulated frame was then analyzed using FOCAS 
and a catalogue of objects produced using the same procedure implemented for the actual 
dataset. Using the parameters corresponding to Fig. 2a (5 < 6 < 34 arcsec, 20 < r& < 23, 
23 < r s < 24), P(f>(<f>) for the simulated faint galaxy images was computed. Again, neither 
the KS nor the \ 2 test rejects the uniform distribution for the artificial data frame and we 
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conclude that our observed non-uniform P(f>(<f>) is not a result of errors in the analysis. 

From the results of our various tests we therefore conclude that the signal observed in 
Fig. 2 is real and not an artifact of the dataset. In the following section we develop a model 
to recreate the observed non-uniform P(f>(<f>) using galaxy-galaxy gravitational lensing. 

3. SIMULATIONS OF GALAXY-INDUCED IMAGE DISTORTIONS 

In order to compare the magnitude of our apparent gravitational lensing signal with 
theoretical expectation, we have carried out a variety of analytic and Monte Carlo sim- 
ulations of weak lensing of galaxy images. To do this we adopt a simple, fiducial galaxy 
model, compute the expected signal using analytic and Monte Carlo simulations, then 
explore the sensitivity of the predictions to changes in our assumptions. In so doing, we 
are implicitly assuming that the majority of our lens galaxies form a one parameter family 
as far as their gravitational properties outside ~ 30 kpc are concerned. Some cosmogonic 
theories posit that the outer parts of ellipticals are similar to the halos of spirals (eg. de 
Zeeuw & Franx 1991), a view that is supported by X-ray observations (eg. Fabbiano 1989). 
We shall therefore not distinguish the minority of ellipticals from the spirals present in our 
lens sample. 

3.1 Galaxy mass model 

In devising an appropriate galaxy model, two simplifications present themselves. 
Firstly, we are quite unconcerned with the details of the galaxy cores as we expressly 
exclude 8^5 arcsec. We can therefore adopt a model galaxy potential that is singular 
as r — > 0. Secondly, almost all of our images are weakly distorted as they lie well outside 
the tangential critical curves formed by the galaxies. This implies that if the ellipticity 
in the potential is e g {r) at radius r, the average image polarisation (p) (r, e g ) will differ 
from that of a circular lens with the same mass contained within a circle of radius r by 
an amount 0(e 2 (p) (r, 0)). As e ^ 0.3 for the bulk of the images, this correction is small 
compared with the uncertainty implicit in the model. We therefore treat the galaxy po- 
tential as circularly symmetric. It should be emphasized at this stage that it is not at 
all likely that the outer parts of galaxies will have had time to relax completely into the 
ellipsoidal shapes associated with the luminous inner regions. In particular, the largest 
equipotentials associated with individual galaxies may be quite non-circular. However, as 
weak lensing depends linearly upon the mass distribution and we are only attempting a 
statistical measurement, our model corresponds to a circularly- averaged mass distribution. 
Note that when an image is distorted by two or more distinguishable lens galaxies, their 
linear effects can be superposed. 

After some experimentation, we have found that a simple model mass distribution for 
the dark matter halos which contains an outer characteristic scale, s, as a free parameter 
is 

V c 2 s 2 

P( r ) = A ^ 9/ 9 9T (3-1) 

1 K ' ^Qr 2 (r 2 + s 2 ) ' 

where V c is the de-projected circular velocity for r « s and Q is Newton's constant. For 
r >> s, the density declines as r -4 and the polarisation signal becomes small. 
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The surface density associated with this space density is given by 

/OO y2 
^ dxp[(R 2 + x 2 ) 1 / 2 ] = ^L_[i _ (1 + X- 2 )- 1 / 2 ] (3.2) 

where X = R/s. The total mass is finite 



M = ^ (3.3) 

and the mass contained within a radius, r, is given by 

V 2 s 

M(r) = ^Tan _1 (r/5). (3.4) 

y 

We shall need the 2-dimensional potential, relative to the potential at some large but finite 
radius, Rl 

$(i?) = -2 dR'^D- = \(1 + X 2 ) 1 / 2 - X - ln[(l + X 2 ) 1 / 2 + 1]1 (3.5) 

J R R' y v J 

where we drop a constant, logarithmic contribution from the outer reference radius which 
does not contribute to the deflection and the polarisation. In order to demonstrate that 
this potential is physically realisable, we compute an associated distribution function in 
Appendix A and demonstrate that it is stable to small perturbations. 

3.2 Image polarisation 

Using the model potential for an individual lens galaxy given by Eq. (3.1), we now 
compute the deflection angle of light rays and the resultant image polarisation. Let the 
lens galaxy have angular diameter distance Dd and source galaxy have angular diameter 
distances D S} Dds as seen from Earth and the lens, respectively. The deflection of a ray 
from the source is then given by 

(eg. Blandford & Narayan 1992; Schneider, Ehlers & Falco 1993). Equating the image 
polarisation, p(X) } to the induced ellipticity we have 

, N D d X d f a(X)\ N , N 

* x > = --irjx{-ir) = '* G M (3 - 7) 

where a positive polarisation corresponds to a net tangential elongation of the images. The 
coefficient po is given by 

2nV? D d D da AQMDdDds ( . Q , 

P ° ~ sD.c* ~ s*D a c* (3 - 8) 
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Fig. 3: Scaled polarisation variation with angle 6 for a given galaxy. G{X) is defined in 
Eq. (3.9) and X = R/s. 

and the function G{X) is 

G{x) = xm + xiy/i • (3 - 9) 

For X « 1, G(X) ~ 1/X, whereas for X » 1, G(X) ~ 2/X 2 . Note also that G(l) = 
0.88 so that the polarisation p ~ po when R ~ s. G{X) is plotted in Fig. 3. Note 
that Eq. (3.6), (3.7) allow us, in principle, to invert a measurement of p(6) to obtain 
the average galaxy potential. However our data are too sparse and our knowledge of the 
redshift distributions and the luminosity variation of galaxy properties is too incomplete 
to make this practical. 

In order to use our model above to determine the expected image polarisation for 
galaxies in an observational sample, we need to relate the luminosities of galaxies to the 
depths of their potential wells. To do this, we use the observations that the local spiral 
galaxy population can be well-described by a luminosity function scaled to a characteristic 
luminosity L* (eg. Loveday et al. 1992) and that these spirals appear to obey a Tully-Fisher 
relation that relates their luminosities, L, to their circular velocities, V c (eg. Aaronson & 
Mould 1983). With this in mind, we now introduce two scaling laws roughly consistent 
with these observations. 

Firstly, we assume that the circular velocity, V c , scales as the fourth root of the total 
luminosity in a given band, in agreement with the Tully-Fisher relation (and also with the 
Faber-Jackson and fundamental plane laws for ellipticals if we treat V c /y/2 as the central 
velocity dispersion [cf. de Zeeuw & Franx 1991]). Introducing a scaling circular velocity, 
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V*, we have 

where the subscript r refers to the r-band and we assume that the emitted spectrum does 
not change with cosmological epoch. A fiducial estimate of V* is 220 km s _1 (eg. Fich & 
Tremaine 1991). 

Our second scaling relation is more of a hypothesis. We assume that the total mass 
to light ratio of a galaxy is a constant independent of its luminosity. In other words we 
suppose that there is a one parameter family of galaxy potentials whose outer radii, s, 
scale as s oc M 1 / 2 oc l}J 2 oc V c 2 . Introducing a scaling radius, s* , we have 

s ( L \ 1//2 / M \ 1//2 

? = (4) r = b?0 • (3 - n| 

Imposing our galaxy scaling laws above, then, we conclude that po does not depend 
explicitly upon the galaxy mass and we can therefore write 

Po =P r- (3.12) 

cD s 

where 

p- = 0.10f V ' .Vf-^-V' (3.13) 
V 220km s _ V \100kpc7 ^ ' 

is a reference polarisation. The total mass of an L* galaxy in this model is 

M* = 1.8 x 10 12 ( ( ^— — ) M Q . (3.14) 

V 220km s"V ViOO/i-^pcy ° V ; 

We now compute the expected mean image polarisation from an observational inves- 
tigation such as the one that we have described in §2. For simplicity, in our fiducial model 
we adopt an Einstein-De Sitter universe with angular diameter distances 

D d = %-a d (l - a 1 / 2 ), D s = ^a s (l-a\' 2 \ D ds = ^(af - a\' 2 ) (3.15) 
-no -"0 -"0 

and expansion factor a djS = (1 + ^d, s ) _1 . 

We must also allow for a spectral (or "K") correction, which can be fairly large in the r- 
band since the spectrum is quite steep at wavelengths blueward of this band (i.e. ~ 6500 A) 
in intermediate redshift spirals. As we believe our lenses have redshifts ~ 0.2 — 0.8 (see 
below), we are concerned with spectral shapes in the range ~ 3500 — 5500 A. Comparing 
with colour measurements of faint galaxies, we find that a = —dhiL v / dhiv ~ 3, very 
approximately. For ellipticals, the effective value of the spectral index, a, is closer to 5 and 
so they should be somewhat rarer in an r-selected faint galaxy sample than they are locally. 
We also ignore the large dispersion in the colours of observed galaxies and treat them as 
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a one parameter family. We can then relate the apparent magnitude to the luminosity 
adopting —18.5 + 5 log h as the absolute r magnitude of an L* galaxy. With this we obtain 

= ^Ho£d_y ^ + ^3 + ai0 0.4(23.9-r) (3 ^ 

at a given emitted frequency. 

In order to proceed further, we need redshift distribution functions for the source and 
lens galaxies as a function of apparent magnitude. Lilly (1993) and Tresse et al. (1993) 
have presented I-selected redshift surveys which show that galaxies in a magnitude range 
comparable to that under consideration here (/ ~ 22 or r ~ 23) have a broad redshift 
distribution extending out to z ~ 1. As the results of our simulation are quite sensitive 
to the form of the redshift distribution used, we adopt a parameterised and normalised 
redshift distribution 

, , 0z 2 e-''/">" 

F "- r >- iw,?).-; • (3 ' 171 

For our fiducial model, we adopt f3 = 1.5. The mode is 1.2zo, the median is z m = IAzq 
and the mean is 1.5^o- More generally, we write 

z = k z [z m + z' m (r- 22)}; 20 < r < 24 (3.18) 

where the constant k z = 0.7 for f3 = 1.5, and the derivative of the median redshift with 
respect to r magnitude is z' m = 0.1 fiducially. 

Next, source galaxies with apparent magnitude r s are selected within an annular ring 
of outer radius # max , and inner radius /# ma x- We now average the function G{X) over this 
ring for lenses of fixed redshift z d and magnitudes r d , obtaining 

fn^ max/ / S dXXG(X) 
G(z d ,r d ) = — 

J OdjOm^l s 

2S '(!-/)-(! + S 2 ) 1 ' 2 + (f + S 2 ) 1 ' 2 (3-19) 



(I-/ 2 ) x 

S + (1 + S 2 f/ 2 



+ 25 In 



s + (p + s 2 y/ 2 



where S = s/ D d 9 ma , x . For S << 1, G oc S 2 ; for S >> 1, the relation becomes asymptoti- 
cally, G ~ 1.7S. 

We next integrate over the source redshift for a fixed lens. 



(p) {z d ,r d ,r a ) 



(^p) G(z d , r d )D d ^ dz a F(z a , r a ) . (3.20) 



Finally, we integrate over the redshift distribution of lens galaxies of magnitude r d to 
obtain 

(p)(r d ,r a )= (^p) ' dz d G(z d ,r d )D d (z d )F(z d ,r d ) J~ ' dz a F(z a ,r a ) (J^j . (3.21) 
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Fig. 4: Theoretical variation of polarisation as a function of source magnitude r s and lens 
galaxy magnitude r& according to Eq. (3.21). 

Note that the polarisation is directly proportional to the assumed velocity dispersion of 
an L* galaxy but depends non-linearly on a single galaxy structural parameter, s* . 

3.3 Fiducial lens model 

To summarise, our fiducial lens model in an Einstein-De Sitter universe uses: 



v* 


= 220kms _1 


* 

s 


= 100/i _1 kp 


a 


= 3 


max 


= 34 arcsec 


/ 


= 0.15 




= 1.5 




= 0.7 




= 0.47 


1 

y 

TO 


= 0.1 



In this model, the median redshift for an r = 22 galaxy is z = 0.47, in rough agreement 
with observations. At this redshift # max corresponds to 118/i -1 kpc. Using this model, 
then, we have evaluated Eq. (3.21) and show the results as a contour plot of (p) as a 
function of r^, r s in Fig. 4. 
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The image polarisation measured from an observational sample is, of course, an av- 
erage over a distribution of source magnitudes, r s , and lens magnitudes, r&. Assuming 
the number counts of galaxies as a function of r-magnitude follow a relation of the form 
= A 10°' 33r , we find the mean image polarisation based upon our analytic model to 
be (p) = 0.011 for lenses and sources corresponding to our principal observational sample 
(20 < r d < 23 and 23 < r s < 24). 

3.4 Calculation of Polarisation 

Using linear theory, we can predict the distribution in image polarisation. As in §2.1 
we write the complex orientation of the image of a source galaxy as x = ee 2 * 1 ^. If the 
polarisation is p, then to linear order the intrinsic source orientation is xo = X ~ P an d the 
effect of the lens galaxy is a simple translation on the x plane. Let this translation be in 
the x direction so that the normalised, observed distribution of orientations is 

fx(Xx,Xy) = f X o(Xx ~P,Xy) (3-23) 

where f x o is the intrinsic distribution in xo- An unbiased estimator of the polarisation is 
then given using 



(x) = / dXxdXyf x (Xx,Xy)Xx 

dxoxdxoyf x o(xo x i XOy )xox +P j dxxdx y fx(Xx,Xy) ( 3 -24) 
= P j dxxdx y fx(Xx,Xy) 
where we have assumed that f x o is isotropic. It is cleanest to estimate p directly using 

p = JdXxdxyf x (Xx,Xy)Xx ^ 3 

JdXxdxyfxi Xx,Xy ) 

However, what we actually measure is the distribution in P(f>(<f>) integrated over ellip- 
ticity. We therefore write 



/[• ir 
de ef x o +pcos2</> / de e— — 

= (£) [l-(p)cos2^( e - 1 >] 



(3.26) 



where the best fit ellipticity distribution, Eq. (2.1), yields 1 'j = 8.0 and (p), the average 
polarisation in the sample, should vary with the lens and source galaxy selection criteria. 

3.5 Monte Carlo Simulations 

To determine the best-fit halo parameters for our lensing model we have constructed 
Monte Carlo simulations of gravitational lensing of background galaxies by foreground 
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galaxies. For every observed object with 20 < r < 24, we assign a galaxy to a random 
location on the Monte Carlo CCD frame (the angular clustering of the observed objects 
is weak, eg. Brainerd, Smail & Mould (1995), and has a negligible effect on the predicted 
signal). The galaxies are then assigned r-magnitudes and ellipticities drawn at random 
from the corresponding observed distributions and position angles uniformly distributed on 
[— 7r/2, 7r/2]. Based on its assigned r-magnitude, each galaxy is then given a redshift chosen 
from the parameterised redshift distribution given by Eq. (3.17), where the parameters 
assumed for the redshift distribution are those summarised by Eq. (3.22). 

The mass distribution of the lens galaxies is modeled according to Eq. (3.1) and we as- 
sume all lenses follow same scaling relations with V* and s* (i.e. Eq. (3.10) and Eq. (3.11)), 
which are the parameters we wish to investigate. For each source galaxy we compute the 
net image polarisation induced by all foreground lens galaxies. Since we are dealing with 
the weak lensing limit, it is sufficient to compute the individual image polarisations due to 
distinguishable lenses and superpose them to obtain the net polarisation. For the adopted 
galaxy redshift distributions we find that roughly a third of the sources are lensed by only 
a single foreground galaxy, another third are affected by two lenses and the remaining third 
encounter 3 or more significant deflectors. 

Having computed the net polarisation of each of the faint galaxy images, we then 
compute P(f>(<f>) for the Monte Carlo images in exactly the same manner as the actual data 
shown in Fig. 2a. By performing many (~ 1000) Monte Carlo realisations for a given pair of 
scaling parameters (V* } s*), we determine a good estimate of the faint galaxy polarisation 
for a particular lensing model. By comparing the predicted and observed (p) and using a 
X 2 minimisation, we can then obtain best-fit parameters for the dark halos of the lenses. 
Results of this procedure are summarised in §4. 

3.6 Sensitivity to model parameters 

As mentioned above, the predicted image polarisation is simply proportional to the 
square of the assumed fiducial circular velocity V*. However, the dependence upon 
the outer radial scale, s* , is more subtle. If we increase s* from its fiducial value to 
~ 300/i -1 kpc, keeping all other parameters in the standard model the same, the polari- 
sation remains constant at 0.011. Only if we reduce s* substantially is there a substantial 
reduction in the measured polarisation. For example, changing s* to 30/i -1 kpc reduces 
(p) from 0.011 to 0.009. The reason for this behaviour is that, within our standard model, 
most of the signal is contributed by lenses that are sufficiently close to the source on the 
sky that the line of sight passes through the isothermal part of the dark halo and, as 
p* oc s _1 for a given velocity dispersion and G oc s, for large halos the average polarisation 
is approximately independent of s. The overall mass to light ratio will, however, increase 
in proportion to s. 

Next we vary the spectral index, a, independently of the other parameters. Changing 
a from 3 to 2 reduces (p) from 0.011 to 0.009. This comparatively weak variation indicates 
that our simplistic model for the K-correction is not a serious concern. 

We now determine the sensitivity of the model predictions to the lens and source 
redshift distributions. Let us keep the median redshift (and its z variation) constant but 
change the shape of the distribution. Reducing the high z tail by increasing f3 from 1.5 
to 2 (a Gaussian) results in a reduction of the mean polarisation from 0.011 to 0.008. 
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Alternatively, we can keep the shape of the redshift distribution unchanged and increase 
the median redshift. In this case, we find that shifting z m from 0.47 to 0.7 results in a 
mean polarisation of 0.015. It turns out that our choice of variation of median redshift 
for the source galaxy distribution with magnitude more or less maximises the predicted 
polarisation and there is only a small sensitivity to z' m . 

Similarly, there is a little sensitivity to the world model. For example, changing from 
an Einstein-De Sitter universe to an empty universe (fio = 0, Dd s = c(2H( i a c i)~ 1 (2 + a s + 
0'd){a'd — O's)) results in a mean polarisation of 0.012. 

Our measurement of the mean image polarisation used a fixed aperture around each 
lens galaxy of 34 arcsec. It turns out that this, too, roughly maximises the signal to 
noise in our simulations. To understand this, first observe that the induced polarisation 
is linear in the mass distribution and so it ought not to matter if several lens galaxies are 
contributing as long as they can be treated as randomly oriented. If we increase # max , then 
the number of lens-source pairs will increase as ^ ax , and the random error will diminish 
as #max- The polarisation signal from these extra background galaxies will diminish as 
^max as l° n g as we are s tiH looking through the dark halos, with # max ^ #c and thus the 
signal to noise ought to improve logarithmically with increasing # max . However, when 
the impact parameters exceed the lens outer radii, s, the signal to noise will deteriorate. 
Consequently, there is little gain in signal to noise possible from increasing the aperture 
size beyond s. For example, increasing # max from 34 to 60 arcsec, leads to a reduction 
in the mean polarisation signal (p) from 0.011 to 0.006, a reduction by a factor 1.8, even 
though the number of galaxies pairs increases by a factor 3.1, which should reduce the 
noise by a factor 0.57, leaving the signal to noise ratio constant. Conversely, reducing # max 
to 15 arcsec, increases (p) by a factor of 2.2, but reduces the number of pairs by 0.18 so 
that the signal to noise diminishes by 0.8. 

4. RESULTS 

Using the functional form for P(f>(<f>) given by Eq. (3.26)and fitting to the observations 
we obtain (p) = 0.011 ± 0.006 (95% confidence bounds) for our fiducial galaxy sample 
(20.0 <r d < 23.0, 23.0 <r s < 24.0) and 5 < 6 < 34 arcsec. As discussed in §2, we cannot 
create such a signal from systematic effects in either our dataset or our analysis. We 
must therefore accept that the signal is real and search for an astrophysical origin. In this 
regard, two processes present themselves as obvious contenders to produce a preferential 
alignment of faint galaxy images around brighter systems: tidally distorted companions or 
gravitational lensing. 

We now show that companion galaxies are unlikely to be the source of the observed 
signal. Limits can be set on the magnitude of the contamination due to tidally-induced 
tangential elongations of genuine satellite galaxies of the lens galaxies (cf. Phillips 1985, 
Tyson 1985) from the clustering strength of the galaxies. Fortunately, the clustering statis- 
tics of the galaxies in this field have already been carefully determined (Brainerd, Smail 
& Mould 1995). We estimate the angular cross correlation function for galaxies in our 
primary data set with 20 < r& < 23, 23 < r s < 24 to be 




(4.1) 



17 



where we have computed the integral constraint directly for our frame and allowed a 15% 
stellar contamination of the object catalogue (see Brainerd, Smail & Mould 1995). Let 
us suppose that that the mean polarisation associated with intrinsic tidal effects is p t . 
Our cursory examination of images of nearby galaxies reveals that \p t \ < 0.1 for satellites 
within a magnitude fainter of their neighbour and separations ~ 100 kpc. The contribution 
of the tidal elongation to the observed polarisation signal should therefore be of order 
1.6w(8 iaiix )pt } which is smaller in modulus than 0.004 for # max = 34 arcsec and of order 
2.3cr below our measured image polarisation. We therefore conclude that a contribution to 
the observed polarisation due to a population of tidally distorted dwarfs can be ruled out 
with moderate confidence. A contrary calculation that illustrates the effect of gravitational 
lensing on the measured autocorrelation function of galaxies at faint magnitudes is given 
in Appendix B. 

We conclude that the most probable cause of the signal reported in §2 is the distortion 
of the distant galaxy images due to gravitational lensing by foreground galaxies. Our 
fiducial model, Eq. (3.22), yields a prediction of (p) = 0.011 for the analytic model and 
(p) = 0.009 ± 0.003 for the Monte Carlo simulations. The theoretical predictions of (p) 
are, thus, in good agreement with each other and with the observational measurement. 
We conclude that we are able to recreate easily the observed lensing signal with a simple, 
physically motivated model, strengthening the case for galaxy-galaxy lensing as the cause 
of the observed image polarisation. 

To determine the level at which we can constrain the parameters of our lensing model, 
we have investigated the radial dependence of the polarisation signal. The image polarisa- 
tion will diminish for lines of sight passing outside the individual dark matter halos of the 
lenses and, so, the scale at which the observed (p) approaches zero is an indication of the 
angular extent of the halos. In Fig. 5a we show (p) for our fiducial subsample of galaxies 
calculated in annuli of inner radius 5 arcsec and outer radius # max , where 15 < # max < 145 
arcsec. These measurements of (p) are, of course, not independent but they do serve to 
illustrate that (p) is significantly non-zero up to scales of 6 max ~ 90 arcsec, suggesting the 
dark halos of the lenses extend to large radii (<; 100/i -1 kpc). 

A more ambitious approach is possible if we attempt to measure the sizes of the 
lens galaxies by determining the variation of polarisation with the differential source-lens 
separation, 6. Again using our fiducial subsample of galaxies, we have computed (p) as a 
function of 6 and the results are exhibited in Fig. 5b. In addition, we have computed the 
mean polarisation of the faint galaxy images in the Monte Carlo simulations as a function 
of 6 for our fiducial L* galaxy model assuming different values of the outer radial scale, 
s* . Results for V* = 220 km s _1 and s* = 20h~ 1 kpc, 100/i _1 kpc, 250h~ 1 kpc are shown 
in Fig. 5b. 

Keeping V* constant at 220 km s _1 , s* was increased incrementally from 10/i -1 kpc to 
250h~ 1 kpc in the Monte Carlo simulations and the variation of (p) with 6 was determined 
for each value of s* . A \ 2 statistic comparing the model predictions to the observed 
(p) in Fig. 5b was then computed. As s* was increased from 10/i -1 kpc, \ 2 decreased 
monotonically from ~ 7.5 and reached a constant value of ~ 3.0 for s* ^> 100/i -1 kpc. 
That is, for values of s* ^> 100/i -1 kpc, the expected variation of (p) is essentially constant, 
as discussed in §3.6. We are, therefore, unable to determine a unique best-fit value of s* , 
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Fig. 5: Angular variation of image polarisation for foreground galaxies with 20 < r& < 
23 and background galaxies with 23 < r s < 24. (a) Variation of (p) with increasing 
annulus outer radius, # max - (b) Variation of (p) with differential lens-source separation, 
6. Theoretical estimates of (p) for fiducial L* galaxy gravitational lenses (see §3) with 
different scaling radii, s* , are also shown. 

but the x 2 minimisation suggests that the halos of the lenses would be at least 100/i -1 kpc 
in radial extent. 

It is clear from Fig. 5b that our fiducial model with s* <; 100/i -1 kpc yields a good fit 
to the observed variation of (p) with source-lens separation. However, it is also interesting 
to investigate the amount by which we can change V* and still obtain a reasonable (p) 
compared to the observations. Keeping s* constant at 100/i -1 kpc, V* was varied in the 
Monte Carlo simulations and, as above, a \ 2 statistic was used to compare the observed 
variation of (p) with 6 to the model predictions. From this we find the upper and lower 
90% confidence bounds on V* to be ~ 304 km s _1 and ~ 116 km s _1 . 
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Finally, it is clear from Eq. (3.17) that the redshift distributions of our fiducial sub- 
samples of lenses (20 < r& < 23) and sources (23 < r s < 24) will have some overlap. We 
expect the measured mean polarisation to be greatest for the most distant sources and 
for the case that all the sources are at higher redshifts than the candidate lenses. We are 
therefore driven to look at using additional information to better separate the foreground 
lenses from the background sources. One obvious possibility is to use the colours of the 
faint sources from our matched g and r-band data. We thus split the source sample with 
23 < r s < 24 into a "red" half [{g-r) > 0.53] and a "blue" half [(g-r) < 0.53] on the basis 
of their (g — r) colours, and compute (p) for the 2 subsamples. Again using 5 < 6 < 34 
arcsec, we find (p) blue = 0.016 ± 0.008 and (p) red = 0.008 ± 0.008 (95% confidence limits). 
In addition, we compute (p) for the red and blue subsamples as in Fig. 5a and in Fig. 6 we 
show the measured variation of (p) with # max . From this figure (p) red appears consistent 
with zero over all scales, while (p) blue is significantly non-zero for # max ^ 60 arcsec and 
there is evidence of a monotonic decrease to zero with increasing # max . While not highly 
significant, these data suggest a higher mean polarisation of the blue images than the 
red images, a result which makes sense if a larger proportion of distant sources are blue, 
star-forming systems. 

5. DISCUSSION 

In this paper we have attempted to measure the induced polarisation of images of 
distant galaxies due to weak gravitational lensing by more nearby galaxies. We have a 
significant detection of this polarisation of (p) = 0.011 ± 0.006 (95% confidence bounds). 
We cannot explain this signal through systematic effects within our dataset and thus 
we believe it is real. In addition, we have presented a fiducial model which is capable 
of reproducing the observed gravitational lensing signal though both analytic and Monte 
Carlo simulations. The simulation results are quite robust to most variations apart from the 
scaling circular velocity, V* , and the details of the redshift distribution. With somewhat 
lower significance we can claim to have measured the decrease in the signal with increasing 
galaxy separation, consistent with a typical galaxy halo size <; 100/i -1 kpc. 

The best-fitting model parameters from the Monte Carlo simulations can be used to 
estimate the masses of the lens galaxies contained within a radius, r. For V* = 220km s _1 
and s* = 100fc -1 kpc, we find M(100/i" 1 kpc) = 8.8 x IQ^H^Mq. For the 90% confidence 
limits on V* derived above this becomes M(100/i _1 kpc) = (S.Slgj) X lO 11 ^ -1 ]^©. Let- 
ting s* — s- oo, we obtain a maximum contained mass of M max (100/i _1 kpc) = (l.l^J'g) X 
10 12 h~ x Mq for our allowed range of V*. From our model calculations, then, we estimate 
an allowed ranee for the masses of the lens halos to be: M ~ l.OtJ^ x IO^H^Mq. 

The typical luminosities of the lens galaxies in our sample, given our fiducial median 
redshift, are Ly ~ 5 X 10 9 h~ 2 Lq. We thus estimate rest-frame mass to light ratios inside 
a radius of 100/i -1 kpc of our lensing galaxies to to be M/L v = 200± 2 1 2 4 ° h(M / L v ) e for 
Oo = 1- A value of M / Ly = 1400/i(M / Lv)q is required for closure density and we, 
therefore, estimate the fraction of the closure density contained in the central regions of 
galaxies to be Q g = 0.14^q'Jq. Our observations thus provide the first evidence from 
weak gravitational lensing of large-scale dark halos associated with individual galaxies. 
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(a) "blue" source galaxies, 23 < r s < 24 
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(b) "red" source galaxies, 23 < r s < 24 
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Fig. 6: Angular variation of image polarisation for the "blue", (g — r) < 0.53, and "red", 
(g — r) > 0.53, source subsamples as a function of outer annular radius, # max , where 
23 < r s < 24. 

Few measurements exist for galaxy masses on the scales probed here. Perhaps the best 
existing estimates of the extent and masses of galaxy halos come from statistical studies of 
satellite galaxies. Zaritsky & White (1994) have analysed a sample of companion galaxies 
to isolated local spirals and find that the typical masses out to 150/i -1 kpc are in the range 
1-2 x 10 12 fc -1 M©, consistent with our findings, M(150/i" 1 kpc) = 1A±\-* x 10 12 fc -1 M©. 
We stress that the assumptions underlying the two techniques are very different which 
makes the close agreement from the two methods very encouraging. 

We conclude that weak gravitational lensing is a viable and potentially powerful probe 
of the outer parts of normal galaxies - regions that are inaccessible to strong lensing 
studies (Breimer & Sanders 1993). That our polarisation measurement has been achieved 
using only a single CCD frame under conditions of modest seeing augurs well for future 
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investigations. When planning such observation a compromise has to be reached between 
area coverage and depth when accruing the source and lens samples. Galaxy counts increase 
at a rat e ~ 2 per magnitude, hence the number of galaxy-lens pairs should quadruple 
for each magnitude and, theoretically, the random error in the polarisation measurement 
should halve, making it advantageous to go fainter. In practice, the limiting magnitude is 
set by our ability to assign accurate shapes to individual source galaxies and we have seen 
that the empirical limit for our data is r ~ 24 (Fig. 2). Note that, in the absence of scale 
evolution, beyond a redshift z ~ 0.5 the image sizes of galaxies are relatively fixed and to 
increase the depth of the sample by ~ one magnitude takes ~ 5 times as much integration. 
Consequently, even when the galaxy images can be measured accurately, galaxy-lens pairs 
are accumulated at about the same rate by taking additional CCD frames as by increasing 
the depth of exposure of an individual image. 

It is clear from our observations and simulations that atmospheric seeing seriously 
degrades the polarisation at faint magnitudes where the signal to noise would otherwise 
increase. This should not be a problem for deep HST images and it is important to repeat 
this measurement using the deepest WFPC2 fields as they become publicly available. 
Using the archival data we anticipate being able to measure image ellipticities accurately 
to r = 26. In this case, the lensing signal should go up by a factor of 3, while the noise 
should go down by a factor of 4 in a given area. Currently, on the order of 10 fields should 
be appropriate for our analysis, so that the signal to noise may improve by a factor of 30 
(although this will depend upon the galaxy redshift distribution). It should, therefore, be 
possible to improve the accuracy of the measurement substantially and constrain both the 
redshift distribution of faint galaxies and the sizes of their halos. 

The eventual limitations to this technique are probably connected with our incomplete 
knowledge of the galaxy redshift distribution and the unknown contribution to the statis- 
tical signal of a minority of unusual lenses and tidal interactions. Despite these concerns, 
the prospect of studying galaxies on a scale where their dynamical times are comparable 
to their ages is exciting. In this regime we can observe the infall of the outer parts of the 
galaxies and thus, in some sense, see the galaxies assembling. This encourages us to devote 
more observational effort to measuring galaxy-induced polarisation. 
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APPENDIX A: Particle distribution function associated with model potential 

In Sec. 3.1, we introduced a simple galaxy model with a density distribution given by 
Eq. (3.1). Associated with this density distribution is the 3D potential 

, f.,GM(r') T/2 

Hr) = -J r dr'—^ = -V 2 

where x = r/s. In order to demonstrate that this potential is physically realisable, we 
compute the associated isotropic distribution function for dark matter particles assuming 
that they are all of the same mass m and ignoring the luminous stars, which ought only 
to contribute at small radius. We denote this distribution function by f(E), where E = 
v 2 /2 + <f) is the specific energy. If we now regard p as a function of <f>, then we can write 

,o 

p(<f>) = 42 1/2 vrm / dE(E - <f>) 1/2 f(E), (A2) 
J<t> 

eg. Binney & Tremaine (1987) This integral equation is easily cast in Abel form and can 
be solved to give 



f(E) = 




where we have imposed the boundary condition that p and its derivatives — > as <f> — > 0. 
This distribution function can be evaluated numerically and is found to decline monoton- 
ically with E. In addition, d 3 p/d<f> 3 < 0. These properties suffice to ensure stability to 
small perturbations (eg. Binney & Tremaine 1987). 



t an x 1 , , o , 

+ -ln(l +x 2 ) 

x A 



(Al) 
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APPENDIX B: Influence of weak lensing on galaxy autocorrelation function 



Curiously, weak gravitational lensing produces an observable effect on the autocor- 
relation function of distant galaxies. It might be thought that the magnification of more 
distant galaxies by intervening lenses would result in a positive contribution to the mea- 
sured autocorrelation function. However, the deflection of the galaxy images also expands 
the separation between galaxies and this turns out to be the larger effect. If the slope of 
the galaxy counts can be expressed as 



where q ~ 0.3, then the surface density of background galaxies at a distance 6 from a given 
lens galaxy will be given be enhanced by a factor 1 + w'(6) where 



and fJ.(6) is the magnification (eg. Narayan 1989). The contribution of gravitational lensing 
to the correlation function turns out to be negative, reducing the very small scale clus- 
tering since the distance between galaxies has been increased. We can relate w'(6) to the 
polarisation, p(6), using Eq. (3. 2), (3. 7), and we find 



This increases from ~ — 1/4 for X << 1 to ~ —1/16X for X >> 1. This ratio is always 
small and so cannot affect the measured image polarisation. However, it can have an 
effect upon the autocorrelation function of galaxies fainter than those scrutinised here by 
canceling a significant fraction of the true correlation on scales <J 10 arcsec. On scales 
of ~ 15-20 arcsec, however, the effect is very small and it is negligible on scales <; 50 
arcsec. Therefore, weak gravitational lensing cannot account for the observed low clustering 
amplitude of very faint galaxies on scales <; 15 arcsec (eg. Brainerd, Smail & Mould, 1995; 
Efstathiou et al. 1991). 



w'{6) ~ (2.5? - 1>(0) 



(B2) 
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Table 1 



faint image 
magnitudes 



number of number of faint 
faint images bright pairs 



X 2 rejection 
confidence level 



KS rejection 
confidence level 



23.0 < r < 24.0 
23.0 < r < 25.0 
23.0 < r < 26.0 



506 

1755 

4303 



3202 

10870 

26412 



98.6% 
97.3% 
19.4% 



99.9% 
99.2% 
60.0% 



FIGURE CAPTIONS 

Fig. 1. Orientation of faint galaxies relative to bright galaxies. 

Fig. 2. Probability distribution P(f>(<f>) of orientation of faint galaxies relative to the direc- 
tions of bright galaxies with projected separations 5 < 6 < 34 arcsec. The bright galaxies 
have 20 < r < 23 and the faint galaxies have magnitudes in the ranges indicated in the 
text and figure panels. For the best case, (a), of a non-uniform P^,(</>), the best fitting 
theoretical cos 2<f> variation is also shown. 

Fig. 3. Scaled polarisation variation with angle 6 for a given galaxy. G{X) is defined in 
Eq. (3.9) and X = R/s. 

Fig. 4. Theoretical variation of polarisation as a function of source magnitude r s and lens 
galaxy magnitude r& according to Eq. (3.21). 

Fig. 5. Angular variation of image polarisation for foreground galaxies with 20 < r& < 
23 and background galaxies with 23 < r s < 24. (a) Variation of (p) with increasing 
annulus outer radius, # max - (b) Variation of (p) with differential lens-source separation, 
6. Theoretical estimates of (p) for fiducial L* galaxy gravitational lenses (see §3) with 
different scaling radii, s* , are also shown. 

Fig. 6. Angular variation of image polarisation for the "blue", (g — r) < 0.53, and "red", 
(g — r) > 0.53, source subsamples as a function of outer annular radius, # max , where 



23 < r s < 24. 
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